Patterns recovered in phylogenomic analysis of Candida auris and close relatives implicate broad environmental flexibility in Candida/Clavispora clade yeasts.

Fungal pathogens commonly originate from benign or non-pathogenic strains living in the natural environment. The recently emerged human pathogen, Candida auris, is one example of a fungus believed to have originated in the environment and recently transitioned into a clinical setting. To date, however, there is limited evidence about the origins of this species in the natural environment and when it began associating with humans. One approach to overcome this gap is to reconstruct phylogenetic relationships between (1) strains isolated from clinical and non-clinical environments and (2) between species known to cause disease in humans and benign environmental saprobes. C. auris belongs to the Candida/Clavispora clade, a diverse group of 45 yeast species including human pathogens and environmental saprobes. We present a phylogenomic analysis of the Candida/Clavispora clade aimed at understanding the ecological breadth and evolutionary relationships between an expanded sample of environmentally and clinically isolated yeasts. To build a robust framework for investigating these relationships, we developed a whole-genome sequence dataset of 108 isolates representing 18 species, including four newly sequenced species and 18 environmentally isolated strains. Our phylogeny, based on 619 orthologous genes, shows environmentally isolated species and strains interspersed with clinically isolated counterparts, suggesting that there have been many transitions between humans and the natural environment in this clade. Our findings highlight the breadth of environments these yeasts inhabit and imply that many clinically isolated yeasts in this clade could just as easily live outside the human body in diverse natural environments and vice versa.

Heritable differences in abundance of bacterial rhizosphere taxa are correlated with fungal necrotrophic pathogen resistance.

Host-microbe interactions are increasingly recognized as important drivers of organismal health, growth, longevity and community-scale ecological processes. However, less is known about how genetic variation affects hosts' associated microbiomes and downstream phenotypes. We demonstrate that sunflower (Helianthus annuus) harbours substantial, heritable variation in microbial communities under field conditions. We show that microbial communities co-vary with heritable variation in resistance to root infection caused by the necrotrophic pathogen Sclerotinia sclerotiorum and that plants grown in autoclaved soil showed almost complete elimination of pathogen resistance. Association mapping suggests at least 59 genetic locations with effects on both microbial relative abundance and Sclerotinia resistance. Although the genetic architecture appears quantitative, we have elucidated previously unexplained genetic variation for resistance to this pathogen. We identify new targets for plant breeding and demonstrate the potential for heritable microbial associations to play important roles in defence in natural and human-altered environments.

A review of the taxonomic diversity, host-parasite interactions, and experimental research on chytrids that parasitize diatoms.

Diatoms (Bacillariophyta) are a major source of primary production on Earth, generating between 1/4 to 1/2 of all oxygen. They are found in almost all bodies of water, the ice of mountains, the arctic and the antarctic, and soils. Diatoms are also a major source of food in aquatic systems, a key component of the silica cycle, and are carbon capturers in oceans. Recently, diatoms have been examined as sources of biofuels, food, and other economic boons. Chytrids are members of the Kingdom fungi comprising, at a minimum, Chytridiomycota, Blastocladiomycota, and Neocallimastigales. Most chytrids are saprobes, plant pathogens, or parasites, and play an important role in aquatic ecosystems. Chytrid parasitism of diatoms has been reported to cause epidemics of over 90% fatality, though most of the information regarding these epidemics is limited to interactions between just a few hosts and parasites. Given the ubiquity of diatoms, their importance in natural and economic systems, and the massive impact epidemics can have on populations, the relative lack of knowledge regarding parasitism by chytrids is alarming. Here we present a list of the firsthand accounts of diatoms reported parasitized by chytrids. The list includes 162 named parasitic chytrid-diatom interactions, with 63 unique chytrid taxa from 11 genera, and 74 unique diatom taxa from 28 genera. Prior to this review, no list of all documented diatom-chytrid interactions existed. We also synthesize the currently known methods of infection, defense, and experiments examining diatoms and chytrids, and we document the great need for work examining both a greater breadth of taxonomic diversity of parasites and hosts, and a greater depth of experiments probing their interactions. This resource is intended to serve as a building block for future researchers studying diatom-parasite interactions and global planktonic communities in both fresh and marine systems.

Five new species of Gloeandromyces (Fungi, Laboulbeniales) from tropical American bat flies (Diptera, Streblidae), revealed by morphology and phylogenetic reconstruction.

This paper describes and illustrates five new species of Gloeandromyces (Ascomycota, Laboulbeniales) associated with tropical American bat flies (Diptera, Streblidae). These are Gloeandromyces cusucoensis sp. nov. from Trichobius uniformis in Costa Rica and Honduras, G. diversiformis sp. nov. from Strebla wiedemanni in Costa Rica, G. plesiosaurus sp. nov. from Trichobius yunkeri in Panama, G. pseudodickii sp. nov. from Trichobius longipes in Ecuador and Panama, and G. verbekeniae sp. nov. from Strebla galindoi in Ecuador and Panama. The description of these five species doubles the number of known species in the genus. Morphological characteristics, host association, and a three-locus (18S nuc rDNA, 28S nuc rDNA, TEF1) phylogenetic reconstruction support placement of these taxa in the genus Gloeandromyces. Three of the new species are polymorphic; they have multiple morphotypes that grow in specific positions on the host integument: G. diversiformis f. diversiformis, f. musiformis, and f. vanillicarpiformis; G. plesiosaurus f. asymmetricus and f. plesiosaurus; and G. verbekeniae f. verbekeniae and f. inflexus. Finally, a dichotomous key to all species and morphotypes is presented.

A combined microscopy and single-cell sequencing approach reveals the ecology, morphology, and phylogeny of uncultured lineages of zoosporic fungi.

Environmental DNA analyses of fungal communities typically reveal a much larger diversity than can be ascribed to known species. Much of this hidden diversity lies within undescribed fungal lineages, especially the early diverging fungi (EDF). Although these EDF often represent new lineages even at the phylum level, they have never been cultured, making their morphology and ecology uncertain. One of the methods to characterize these uncultured fungi is a single-cell DNA sequencing approach. In this study, we established a large data set of single-cell sequences of EDF by manually isolating and photographing parasitic fungi on various hosts such as algae, protists, and micro-invertebrates, combined with subsequent long-read sequencing of the ribosomal DNA locus (rDNA). We successfully obtained rDNA sequences of 127 parasitic fungal cells, which clustered into 71 phylogenetic lineages belonging to seven phylum-level clades of EDF: Blastocladiomycota, Chytridiomycota, Aphelidiomycota, Rozellomycota, and three unknown phylum-level clades. Most of our single cells yielded novel sequences distinguished from both described taxa and existing metabarcoding data, indicating an expansive and hidden diversity of parasitic taxa of EDF. We also revealed an unexpected diversity of endobiotic Olpidium-like chytrids and hyper-parasitic lineages. Overall, by combining photographs of parasitic fungi with phylogenetic analyses, we were able to better understand the ecological function and morphology of many of the branches on the fungal tree of life known only from DNA sequences. IMPORTANCE Much of the diversity of microbes from natural habitats, such as soil and freshwater, comprise species and lineages that have never been isolated into pure culture. In part, this stems from a bias of culturing in favor of saprotrophic microbes over the myriad symbiotic ones that include parasitic and mutualistic relationships with other taxa. In the present study, we aimed to shed light on the ecological function and morphology of the many undescribed lineages of aquatic fungi by individually isolating and sequencing molecular barcodes from 127 cells of host-associated fungi using single-cell sequencing. By adding these sequences and their photographs into the fungal tree, we were able to understand the morphology of reproductive and vegetative structures of these novel fungi and to provide a hypothesized ecological function for them. These individual host-fungal cells revealed themselves to be complex environments despite their small size; numerous samples were hyper-parasitized with other zoosporic fungal lineages such as Rozellomycota.

Comparative genomics and phylogenomic investigation of the class Geoglossomycetes provide insights into ecological specialization and the systematics of Pezizomycotina.

Despite their global presence and ubiquity, members of the class Geoglossomycetes (Pezizomycotina, Ascomycota) are understudied systematically and ecologically. These fungi have long been presumed saprobic due to their occurrence in or near leaf litter and soils. Additionally, they lack an apparent association with other organisms, reinforcing this perception. However, observations of sporocarps near ericaceous shrubs have given rise to an alternative hypothesis that members of Geoglossomycetes may form ericoid mycorrhizae or ectomycorrhizae. This claim, however, has yet to be confirmed via microscopy or amplicon-based studies examining root communities. As a result, our current understanding of their ecology is based on cursory observations. This study presents a comparative analysis of genomic signatures related to ecological niche to investigate the hypothesis of an ericoid mycorrhizal or ectomycorrhizal ecology in the class. We compared the carbohydrate-active enzyme (CAZyme) and secondary metabolite contents of six newly sequenced Geoglossomycetes genomes with those of fungi representing specific ecologies across Pezizomycotina. Our analysis reveals CAZyme and secondary metabolite content patterns consistent with ectomycorrhizal (EcM) members of Pezizomycotina. Specifically, we found a reduction in CAZyme-encoding genes and secondary metabolite clusters that suggests a mutualistic ecology. Our work includes the broadest taxon sampling for a phylogenomic study of Pezizomycotina to date. It represents the first functional genomic and genome-scale phylogenetic study of the class Geoglossomycetes and improves the foundational knowledge of the ecology and evolution of these understudied fungi.

Mr. Toad's Wild Fungi: Fungal Isolate Diversity on Colorado Boreal Toads and their Capacity for Pathogen Inhibition.

The amphibian skin pathogen Batrachochytrium dendrobatidis (Bd) has caused an ongoing biodiversity crisis, including in the locally endangered Colorado boreal toad (Anaxyrus boreas boreas). Although researchers have investigated the bacteria living on amphibian skin and how they interact with Bd, there is less information about fungal community members. This study describes (1) the diversity of culturable fungi from boreal toad skin, (2) which subset of these isolates is Bd-inhibitory, and (3) how Bd affects these isolates' growth and morphology. Most isolates were from the orders Capnodiales, Helotiales, and Pleosporales. Of 16 isolates tested for Bd-inhibition, two from the genus Neobulgaria and three from Pseudeurotium inhibited Bd. Fungal growth in co-culture with Bd varied with weak statistical support for Neobulgaria sp. (isolate BTF_36) and cf Psychrophila (isolate BTF_60) (p-values = 0.076 and 0.092, respectively). Fungal morphology remained unchanged in co-culture with Bd, however, these results could be attributed to low replication per isolate. Nonetheless, two fungal isolates' growth may have been affected by Bd, implying that fungal growth changes in Bd co-culture could be a variable worth measuring in the future (with higher replication). These findings add to the sparse but growing literature on amphibian-associated fungi and suggest further study may uncover the relevance of fungi to amphibian health and Bd infection.

Diploid-dominant life cycles characterize the early evolution of Fungi.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Protocol for single-cell isolation and genome amplification of environmental microbial eukaryotes for genomic analysis.

We describe environmental microbial eukaryotes (EMEs) sample collection, single-cell isolation, lysis, and genome amplification, followed by the rDNA amplification and OTU screening for recovery of high-quality species-specific genomes for de novo assembly. These protocols are part of our pipeline that also includes bioinformatic methods. The pipeline and its application on a wide range of phyla of different sample complexity are described in our complementary paper. In addition, this protocol describes optimized lysis, genome amplification, and OTU screening steps of the pipeline. For complete details on the use and execution of this protocol, please refer to Ciobanu et al. (2021).

Improving Fungal Cultivability for Natural Products Discovery.

The pool of fungal secondary metabolites can be extended by activating silent gene clusters of cultured strains or by using sensitive biological assays that detect metabolites missed by analytical methods. Alternatively, or in parallel with the first approach, one can increase the diversity of existing culture collections to improve the access to new natural products. This review focuses on the latter approach of screening previously uncultured fungi for chemodiversity. Both strategies have been practiced since the early days of fungal biodiscovery, yet relatively little has been done to overcome the challenge of cultivability of as-yet-uncultivated fungi. Whereas earlier cultivability studies using media formulations and biological assays to scrutinize fungal growth and associated factors were actively conducted, the application of modern omics methods remains limited to test how to culture the fungal dark matter and recalcitrant groups of described fungi. This review discusses the development of techniques to increase the cultivability of filamentous fungi that include culture media formulations and the utilization of known chemical growth factors, in situ culturing and current synthetic biology approaches that build upon knowledge from sequenced genomes. We list more than 100 growth factors, i.e., molecules, biological or physical factors that have been demonstrated to induce spore germination as well as tens of inducers of mycelial growth. We review culturing conditions that can be successfully manipulated for growth of fungi and visit recent information from omics methods to discuss the metabolic basis of cultivability. Earlier work has demonstrated the power of co-culturing fungi with their host, other microorganisms or their exudates to increase their cultivability. Co-culturing of two or more organisms is also a strategy used today for increasing cultivability. However, fungi possess an increased risk for cross-contaminations between isolates in existing in situ or microfluidics culturing devices. Technological improvements for culturing fungi are discussed in the review. We emphasize that improving the cultivability of fungi remains a relevant strategy in drug discovery and underline the importance of ecological and taxonomic knowledge in culture-dependent drug discovery. Combining traditional and omics techniques such as single cell or metagenome sequencing opens up a new era in the study of growth factors of hundreds of thousands of fungal species with high drug discovery potential.

A single-cell genomics pipeline for environmental microbial eukaryotes.

Single-cell sequencing of environmental microorganisms is an essential component of the microbial ecology toolkit. However, large-scale targeted single-cell sequencing for the whole-genome recovery of uncultivated eukaryotes is lagging. The key challenges are low abundance in environmental communities, large complex genomes, and cell walls that are difficult to break. We describe a pipeline composed of state-of-the art single-cell genomics tools and protocols optimized for poorly studied and uncultivated eukaryotic microorganisms that are found at low abundance. This pipeline consists of seven distinct steps, beginning with sample collection and ending with genome annotation, each equipped with quality review steps to ensure high genome quality at low cost. We tested and evaluated each step on environmental samples and cultures of early-diverging lineages of fungi and Chromista/SAR. We show that genomes produced using this pipeline are almost as good as complete reference genomes for functional and comparative genomics for environmental microbial eukaryotes.

Fungi in the Marine Environment: Open Questions and Unsolved Problems.

Terrestrial fungi play critical roles in nutrient cycling and food webs and can shape macroorganism communities as parasites and mutualists. Although estimates for the number of fungal species on the planet range from 1.5 to over 5 million, likely fewer than 10% of fungi have been identified so far. To date, a relatively small percentage of described species are associated with marine environments, with ∼1,100 species retrieved exclusively from the marine environment. Nevertheless, fungi have been found in nearly every marine habitat explored, from the surface of the ocean to kilometers below ocean sediments. Fungi are hypothesized to contribute to phytoplankton population cycles and the biological carbon pump and are active in the chemistry of marine sediments. Many fungi have been identified as commensals or pathogens of marine animals (e.g., corals and sponges), plants, and algae. Despite their varied roles, remarkably little is known about the diversity of this major branch of eukaryotic life in marine ecosystems or their ecological functions. This perspective emerges from a Marine Fungi Workshop held in May 2018 at the Marine Biological Laboratory in Woods Hole, MA. We present the state of knowledge as well as the multitude of open questions regarding the diversity and function of fungi in the marine biosphere and geochemical cycles.

Leveraging single-cell genomics to expand the fungal tree of life.

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

Harnessing the power of phylogenomics to disentangle the directionality and signatures of interkingdom host jumping in the parasitic fungal genus Tolypocladium.

Host specialization is common among parasitic fungi; however, there are examples when transitions in host specificity between disparately related hosts have occurred. Here, we examine the interkingdom host jump from insect pathogenicity and mycoparasitism in Tolypocladium. Previous phylogenetic inferences made using only a few genes and with poor support reconstructed an ancestral character state of insect pathogenesis, a transition to mycoparasitism, and reversions to insect pathogenesis. To further explore the directionality and genes underlying the transitions in host, we sequenced two additional species of Tolypocladium (T. capitatum and T. paradoxum) and used phylogenomics to compare two insect pathogens and two mycoparasites. Our whole-genome-scale analysis suggests that the diversification of Tolypocladium species happened relatively quickly and that the truffle parasites form a monophyletic, derived lineage within the genus that is the result of a single ecological transition or host jump from insects to fungi. A significant amount of gene tree/species tree discordance occurs within the data set, and we infer this to be the product of both an historical hybridization event and incomplete lineage sorting that was likely because of the rapid diversification of the clade. Furthermore, comparative genomic analyses revealed a set of genes that are exclusive to the mycoparasitic species. These potentially mycoparasitic gene clusters were characterized by a reduced proportion of secreted proteins when compared with entomopathogen-enriched genes and involved the reshaping of the fungal secretome in the ecological context of mycoparasitism.

A phylogenetically-based nomenclature for Cordycipitaceae (Hypocreales).

The ending of dual nomenclatural systems for pleomorphic fungi in 2011 requires the reconciliation of competing names, ideally linked through culture based or molecular methods. The phylogenetic systematics of Hypocreales and its many genera have received extensive study in the last two decades, however resolution of competing names in Cordycipitaceae has not yet been addressed. Here we present a molecular phylogenetic investigation of Cordycipitaceae that enables identification of competing names in this family, and provides the basis upon which these names can be maintained or suppressed. The taxonomy presented here seeks to harmonize competing names by principles of priority, recognition of monophyletic groups, and the practical usage of affected taxa. In total, we propose maintaining nine generic names, Akanthomyces, Ascopolyporus, Beauveria, Cordyceps, Engyodontium, Gibellula, Hyperdermium, Parengyodontium, and Simplicillium and the rejection of eight generic names, Evlachovaea, Granulomanus, Isaria, Lecanicillium, Microhilum, Phytocordyceps, Synsterigmatocystis, and Torrubiella. Two new generic names, Hevansia and Blackwellomyces, and a new species, Beauveria blattidicola, are described. New combinations are also proposed in the genera Akanthomyces, Beauveria, Blackwellomyces, and Hevansia.

The genome of an intranuclear parasite, Paramicrosporidium saccamoebae, reveals alternative adaptations to obligate intracellular parasitism.

Intracellular parasitism often results in gene loss, genome reduction, and dependence upon the host for cellular functioning. Rozellomycota is a clade comprising many such parasites and is related to the diverse, highly reduced, animal parasites, Microsporidia. We sequenced the nuclear and mitochondrial genomes of Paramicrosporidium saccamoebae [Rozellomycota], an intranuclear parasite of amoebae. A canonical fungal mitochondrial genome was recovered from P. saccamoebae that encodes genes necessary for the complete oxidative phosphorylation pathway including Complex I, differentiating it from most endoparasites including its sequenced relatives in Rozellomycota and Microsporidia. Comparative analysis revealed that P. saccamoebae shares more gene content with distantly related Fungi than with its closest relatives, suggesting that genome evolution in Rozellomycota and Microsporidia has been affected by repeated and independent gene losses, possibly as a result of variation in parasitic strategies (e.g. host and subcellular localization) or due to multiple transitions to parasitism.

Morphological, molecular, and ultrastructural characterization of Rozella rhizoclosmatii, a new species in Cryptomycota.

Rozella is a genus of unwalled endoparasites of a variety of hosts including Oomycota (Stramenopiles), Blastocladiomycota and Chytridiomycota (Fungi), and one green alga (Coleochaete, Chlorophyceae). It currently includes more than 20 formally described species, and no new species of Rozella have been described since 1987. We discovered a new Rozella species parasitizing Rhizoclosmatium globosum (Chytridiales, Chytridiomycota) and investigated its morphology, ultrastructure, and phylogenetic position. Herein named as Rozella rhizoclosmatii sp. nov., the organism induces hypertrophy of the host. Its zoospore is ultrastructurally similar to that of Rozella allomycis, although it has a unique zoospore ultrastructural feature, a lattice of perpendicular rods about the nucleus. The 18S rDNA molecular sequence of R. rhizoclosmatii is similar to that of the previously sequenced 'Rozella ex Rhizoclosmatium'. This is the first study to inclusively characterize a new species of Rozella with morphological, ultrastructural and molecular data. As this is only the second Rozella species to be examined ultrastructurally, and because it is parasitic on a member of Chytridiomycota and not Blastocladiomycota, this research supports the conservative nature of zoospore ultrastructure to help define the genus.

Ectomycorrhizal ecology is imprinted in the genome of the dominant symbiotic fungus Cenococcum geophilum.

The most frequently encountered symbiont on tree roots is the ascomycete Cenococcum geophilum, the only mycorrhizal species within the largest fungal class Dothideomycetes, a class known for devastating plant pathogens. Here we show that the symbiotic genomic idiosyncrasies of ectomycorrhizal basidiomycetes are also present in C. geophilum with symbiosis-induced, taxon-specific genes of unknown function and reduced numbers of plant cell wall-degrading enzymes. C. geophilum still holds a significant set of genes in categories known to be involved in pathogenesis and shows an increased genome size due to transposable elements proliferation. Transcript profiling revealed a striking upregulation of membrane transporters, including aquaporin water channels and sugar transporters, and mycorrhiza-induced small secreted proteins (MiSSPs) in ectomycorrhiza compared with free-living mycelium. The frequency with which this symbiont is found on tree roots and its possible role in water and nutrient transport in symbiosis calls for further studies on mechanisms of host and environmental adaptation.

Fine-scale spatial genetic structure of a fungal parasite of coffee scale insects.

The entomopathogenic fungus Lecanicillium lecanii persists in a highly dynamic network of habitat patches (i.e., a metapopulation) formed by its primary host, the green coffee scale Coccus viridis. Lecanicillium lecanii is an important biological control of both C. viridis and the coffee rust, Hemileia vastatrix. Successfully managing this biocontrol agent will depend on an increased understanding of the characteristics of its dispersal, as migration between occupied and unoccupied patches is essential for the persistence of this metapopulation. In the present study, we employ a population genetics approach, and show that in our study system, a coffee farm in the Soconusco region of southern Mexico, L. lecanii is characterized by clear spatial genetic structure among plots within the farm but a lack of apparent structure at smaller scales. This is consistent with dispersal dominated by highly localized transport, such as by insects or rain splash, and less dependence on longer distance dispersal such as wind transport. The study site was dominated by a few multi-locus microsatellite genotypes, and their identities and large-scale locations persist across both study years, suggesting that local epizootics (outbreaks) are initiated each wet season by residual propagules from the previous wet season, and not by long-distance transport of propagules from other sites. The index of association, a measure of linkage disequilibrium, indicates that epizootics are primarily driven by asexual, clonal reproduction, which is consistent with the apparent lack of a teleomorph in the study site and the presence of only a single mating type across the site (MAT-1-2-1). Although the same predominant clonal genotypes were found across years, a drastic difference in genotypic diversity was witnessed across two sites between the two years, suggesting that interclonal selection was occurring. In light of the dispersal limitation of L. lecanii, spatial structure may be an essential axis of management to ensure the persistence of L. lecanii and preserve the ecosystem services provided by this versatile biocontrol agent in this and similar coffee farms.